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ABSTRACT 


Through the Optimum Path Aircraft Routing System (OPARS), 
Fleet Numerical Oceanography Center (FNOC) has committed it- 
self to providing a computerized flight planning service 
remotely accessible via dial-up communications lines. The 
question arises as to whether the proposed number of tele- 
phone lines -will be adequate to provide a level of service 
previously provided by the Lockheed Jetplan system. This 
study provides a detailed analysis of the response delays for 
the OPARS flight plan system. In addition, estimates are 
given of communication requirements when various levels of 
demand prevail, and under conditions in which the FNOC computer 


1s busy or idle. 
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NOTATIONS AND ABBREVIATIONS 


m-Server queueing system where A represents the 
the interarrival distribution and B represents the 
service distribution 

Central Memory 


Central Processor Unit 


Denotes r stage Erlang distribution 


E[T(Setup) ] Expected service time of IRG sub-system 


E{LT(Input queue) ] Expected service time of Input queue sub-system 


E[T(Run) J] 


E[T(Service) ] 


FNOC 


Expected service time of execute sub-system 
Expecter service time for OPARS system 
Fleet Numerical Oceanography Center 


Denotes general distribution 


Arrival rate of OPARS jobs 

eri valerave Cf higher priority jobs 

Sievece rave 

Weneves Expenenvial distribution 

Number of customers in the system at time t 

Long run average number of customers in the system 
Optimum Path Aircraft Routing System 


Pale? LOrmula svacre prebabilities for states 
Set ge ccs. guy ol 


UCPeeZauieon factor 





I. SUMMARY 


In the proposed OPARS flight plan system, users at 
37 remote terminals (Figure 1) connect directly to the 
Fleet Numerical Oceanography Center (FNOC) computer system 
via 11 telephone lines. The total amount of time that a line 
is busy, the total OPARS service time, can be modeled as 
the sum of three sub-system service times: 

1. An interactive program setup time T(setup); 

2. A delay that the OPARS program experiences in 

the FNOC computer's input queue T(input queue) ; 

fee an OPARS program run time T(run). 
The expected values of these individual sub-system service 
times can be computed and then summed to provide a mean 
or expected service time for the entire OPARS request as 
POLL OwW Ss : 

Meeoservice) ] = E[T(setup)] + E[T(input queve)] + E[T(run) ]. 
From this, Erlang's formula equation (1) can be employed to 
compute individual long-run state probabilities (the 


probability that i lines are busy). p;, and ultimately the 


« 
é4 


erobability, of having all eleven lines busy. 


Pia 
imesupperse Of the analytic model, a computer program 

was written to simulate the entire OPARS flight plan request 

Bmecess, {rom initial dial-up through program completion. 

The simulation also includes the arrival and servicing of other 


FNOC computer programs which interfere with the processing 


of the OPARS program. 





meee tl. Remove Terminal Sites 


MOFFET NAS (2) 
ALAMEDA NAS (2) 
JACKSONVILLE NAS (2) 
BRUNSWICK NAS 

NORTH ISLAND NAS (3) 
NORFOLK (3) 

BARBERS POINT NAS (2) 
CHANUTE WY 

NEW ORLEANS 

DETROIT 

SOUTH WEYMOUTH 
WILLOW GROVE 
GLENVIEW 

WHIDBY IS. 


POINT MAGU 


feos, COast Guard Station 


PATUXENT RIVER 
MEMPHIS 

Choa Ta eD 

CHERRY POINT 

EL TORO 

WASHINGTON DC (FAA) 
ieee ly CC .G.) 
KODIAK (C.G.) 
SACRAMENnOnm (CG .G. ) 

ST. PETERSBERG (C.G.) 
Pee OC mace. GC. ) 


MOBILE (C.G.) 


VAS eNe TON EDC (C.G.) 


BARBERS PT (C.G.) 





Both the analytic model and the simulation, computed 
state probabilities for OPARS demand rates of 2, 4, 6, ...20 
per hour and under conditions in which the FNOC computer is 
busy or idle. The results clearly indicate that the ll 
telephone lines can handle demand rates up to three times 
as great as those currently being experienced by the Lockheed 
Jetplan system before Pip the probability of all lines 


busy, exceeds, .0Q5. 
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ie ave ODUCTION 


Fleet Numerical Oceanography Center (FNOC) is currently 
testing and evaluating a computerized flight plan systen, 
referred to, for short, as OPARS. This sytem, developed to 
replace the Lockheed Jetplan flight plan sytem, provides users 
at remote sites with direct access to the FNOC computer 
via ll telephone lines. The purpose of this study is to 
determine if the intended number of telephone lines would 
be adequate to ensure a low probability of having all lines 
busy. 

The number of lines busy at any time t, {N(t)|t>o} can be 
modeled as aBirth-Death process with inter-arrival times 
that are assumed to be independent and exponentially 
Gistributed but with service times that clearly are not. 
Fortunately, in a communications system characterized by 
calls that arrive in a@ stationary or time-homogenous Poisson 
Process of rate Aand vie for n lines with no queue (blocked 
Paes lost), the long run probability of having i lines busy, 


Peaecan be computed from Eriang's formula. 


i 
p. =)(At) Ji! O<i<N a) 
aa aa 
Say! 
j= 
1>N 
0 


This formula has the surprising feature of being valid for 
any service time distribution. Indeed. given the above 


assumption, one need only obtain the mean service time, t, 


| ee 





in order to calculate p,. The Reyer uo, Ol some ClrTory in 

this study was to characterize the OPARS system in such a 

way so that a mean service time could be derived under various 
Seeracving conditions. 

This paper begins with some introductory material 
concerning the FNOC computer systems,the OPARS flight plan 
system and historical usage of the Lockheed Jetplan system. 
Section IV is a detailed description of the analytic model 
meeoomwana 2a discussion of the supporting siumulation program. 
Finally, Section V contains the tabled results and an analysis. 

Notation used in this paper is consistent with that found 
in Queueing Theory literature and will be explained whenever 
Mameoawieced. itn addition, a listing and discription of 


all notation and abbreviations can be found on page seven. 


na 





IIIT. BACKGROUND 


fe PNOC COMPUTER SYSTEM 

PNOC currently has four computers in service at their 
Monterey facility, but only one, eneaee CG as HAL 1s 
accessible to the remote OPARS terminals. HAL is a CDC-6600 
computer with two central processors and 330 K of octal 
central memory. HAL is accessed within FNOC by approximately 
twenty-five interactive terminals as well as by a batch 
job system that routinely processes hundreds of programs a 
day. Since HAL alone is involved in providing service to 
OPARS remote terminalis, we will only consider it in the 
Beeoeussion. 

i. HAL Batch System 

HAL's batch system is composed of a priority ranked 

"first-in, first-out" input queue, which can be considered 
exterior to the computer, and a priority ranked execute 
queue from which jobs are selected by the scheduler for 
processing (Figure 2). When a job enters the system, it 
first enters the input queue, at which the job's user-assigned 
Meormey establishes its initial position. While in the 
queue the job slowly accrues additional "wait time” 
priority which ensures that even the lowest priority jobs 
mee eventually run. The largest number of daily jobs are 
restricted to priorities of 1, 2 or 3 (3 being the highest) 
bao priorities up to 7/7 can be assigned. The priority of 
OPARS jobs is fixed automatically at 60 which immediately Puts 


it ahead of all but a few other Jobs in the queue. 


nn 





When space is available within the computer the scheduler 
removes the first job (with respect to priority then to 
time-of-arrival) from the input queue and moves it to the 
execute queue. Upon entering the execute queue the job 
Goses all of its previously gained "wait-time" priority but 
immediately begins to gain it back as the job waits to be 
processed. In this queue the jobs with the higher assigned 
Memerweres gain this additional priority faster than those 
with low assigned priority. This mechanism results in higher 
priority jobs getting selected for processing sooner than 
low priority jobs. When selected by the scheduling routine, 
the program is moved into central memory and it is processed 
Pommamullt length of time or until it attempts to access 
a device (disk, tape, etc.) which is unavailable. When this 
event occurs the job is removed form central memroy and 
returned to the execute queue, having its accured priority 
Reset to zero. The process continues in this manner until 
the program is completed and it is transferred to the out- 
put queue. 

2. Current FNOC Workload 

From available central processor and central memory 
utilization profiles (Figure 3) it can be readily seen that 
resource demands on HAL can be separated into two states: 

A Busy or High Demand state in which the computer is operating 
at maximum capacity, and an Idle or Low Demand state in which 
most of the resources are immediately available. The High 
Demand state is characterized by a full execute queue, and a 


lengthy input queue, while the Low Demand state exhibits an 
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execute queue with only an occasional job in it and an in- 
put queue which is empty. Fortunately the transition 
meemroamoelween these states is quite brief and therefore 
~ereeaolce lon state is not needed. 

The High Demand period coincides, not surprisingly, 
With normal working hours except that it extends well into 
the evening, often as late as midnight. Any jobs entering 
the system during this period must be delayed for a time 
in the input queue before being run. During the idle period, 
on the other hand, jobs pause only momentarily in the input 


Queue before being run. 


B. OPARS FLIGHT PLAN SYSTEM 

The OPARS flight plan system provides the user with 
an optimum (with respect to time) route of flight given forecast 
weather and user inputs such as aircraft type and. take off 
weight. For flights within areas that have a defined 
routing structure the optimization is fairly easy because 
only a few routes are candidates for the optimum. Over 
feaeer Or point to point flight plans are more difficult 
because there is, in theory, an uncountably infinite number 
ef available routes. 

eaemweware Gwo major sections cf the OPARS program: an 
maeeractive portion and the optimizaticn program. When a 
user at a remote site requests an OPARS flight plan one 
first sets up the program with the Interactive Response 
Generator (IRG). Using a query and response technique, the 


IRG prepares the main program by obtaining required 
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information from the user. During tnis session the IRG 
mommouschecking input information for validity, but only 
mom tormat. For example, an entry of ABCE, as the four- 
iiweer identification code for the destination airfield, 
would be accepted even though no such airfield code 
exists, When all of the required information is entered, 
the IRG allows the user to review and change any input 
Mmemmeecessary. Then, upon command, the program is put into 
the HAL computer batch system. 

As discussed in the previous section, the OPARS job 
enters the input queue along with all other FNOC jobs and 
[eemeooe ibs curn. Hereafter the OPARS program is handled 
memeeiiys OUher batch program with one exception. Upon 
cememet2On the flight plan is returned automatically to the 
terminal where the request was entered. If for some 
Meee) the line has been disconnected, the flight plan is 
placed into an output file from which it can be retrieved 


later. 


C. HISTORICAL USAGE/PROJECTED DEMAND 

In an effort to evaluate expected demand for the OPARS 
flight plan, several methods were employed. They were: 

1) A linear regression model of historical usage by 

the Navy to estimate demand through 1981; 

2) expectations of remote site users; 

3) Lockheed's records of previous use. 
None of these methods alone provided the required level of 
mSeeuracy and, only by combining the three, could a reasonable 


estimate be made. 
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The Navy's own records provided the best overall 
/mieermation. This information, shown in Figure 4, consists 
of monthly totals of Navy requests for the Jetplan. The 
trend is clearly an increase in usage and a simple linear 
regression yielded the projected demands shown in Figure 5. 
Miese monthly totals, unfortunately did not provide any 
information as to how requests varied during the day. 

of — 1978 TCS: 


Month(per day) Month(per day) Month(per day) 
J AN goo0( 26.6 TAG sic ves CO; 






FEB 1012(34.9) Oe >) 1926(68.8) 
MAR 1076(34.7) 1856(59.9) 2107(68) 
APR 1169(38.9) 1319(44) 1887( 62.9) 
fem 1103(35.6) 1490( 48.1) 2119(68.4) 
mt 1174(39.1) 1224(40.8) 2081(69.4) 
te 1 711(55.2) LEGG) 2011(64.9) 
AUG 1066(34.4) 1365(44) ZAG (2 257) 
ear Wi 3253) 1226(40.9) eee ee. 9) 
Sere 1142(36.8) 1415(45.6) 

NOV 1103(36.8) 1321(44) 

wee 1107(35.7) 1386(44.7) 


Figure 4. Historical Usage of Lockheed's Jetplan 


be INTERCEPT Deis 0S JUL 80 ESTIMATE 80.4 /DAY 
SUOPE 1.2 JUL 81 ESTIMATE 94.8 /DAY 


Con. CORRE.  .826 


Figure 5. Linear Regression of Usage Data 
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iemerder tO Obtain the distribution of daily requests, 
expected demand information was requested from all remote 
peeer Users. The resulting information suffered from 
impreeision and a lack of completeness (only about half of 
the remote sites responded). 

As a final effort, an attempt was made to secure actual 
daily records from the Lockheed Corporation. The Lockheed 
personnel were unable to provide me with actual data, but 
over the course of several interviews a picture of the 
daily demand distribution was pieced together. 

The results of these data collection efforts appear in 
Figure 7. These graphs show the projected demand for 
July 1980 and July 1981 distributed using the composite 


results of Lockheed information and user expectations. 
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ive) LAE MODEL 


We start the modeling discussion by reviewing 
some basic assumptions about the system. 

Mmirse and foremost is the assumption concerning 
meme=nomogceneous Poisson arrivals for OPARS requests. 
Without this assumption, the Erlang formula (1) would not 
be valid and a much more complex model would be required. 

secondly, the service times for the three sub-systems 
are assumed to be mutually independent. 

Mmemalliy.s as discussed in the last section, the 
conditions for the computers workload will be divided in- 
BO two distinct periods. These are the High Demand period 
with a nonempty input queue, and the Low Demand period 
where the OPARS requests bypass the input queue and go 


geaeectly into the computer. 


A. THE ANALYTIC MODEL 

Since the main thrust of the analytic effort is to 
obtain an overall mean service time, an immediate 
simplification would be to somehow break the OPARS system 
into two or more sub-systems whose service times would 
be more easily obtained. With this in mind, we start by 
separating the service into the interactive sub-system 
and the batch sub-system. The latter is still rather 
formidable because of the nature of the input and 
execute queues and the complexity of the scheduling 


routine. The final step is to reduce the batch system 





into an input queue sub-system and an execution sub-system. 
The reasons for this may not be readily apparent as there 
are other modeling techniques such as an M|G|m multiserver 
queueing system with priorities, which would handle the 
system as a whole. Such a system would require that each 
job priority have a mean service time. Unfortunately, there 
is no correlation between a jots priority and its core 
and central processor requirements. Additional unrealistic 
assumptions would be required which would, in the end, 
reduce the validity of the result. 
imine interactive Sub-system 

First to be considered is the interactive sub-system. 
An estimate for the mean time for this sub-system is 
readily comput ble from the available, albeit scarce, data. 
This data(Figure 7) was obtained by timing FNOC technicians 
on trial OPARS runs. The mean of these data is 5.1 minutes 
and will be assumed to be the same during both High and 
Low Demand periods. 

2. Input Queue Sub-System 

The next sub-system to model is the input queue sub- 
system. Under low demand conditions the mean service time 
is, by definition, equal to zero. During the High Demand 
Pemrvod if is mot quite so simple, Again, by previous 
definition, the High Demand period has a continuously non- 
empty input queue. This very simple assumption suggests that 
the input queue may be modeled as a single-server queueing 
system by itself, the server being the first position in the 


queue. One must merely postulate the interarrival process 
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Figure 7. Interactive Setup Times 


and the service process, the latter representing the times 
between successive departures from the input queue. If we 
can now model the arrival rate as Poisson and the service 
distribution as exponential we have an M/M/1 queueing 
Eeowem ranked by priority for which there is a simple 
elosed form solution for expected waiting times. The 
interested reader is urged to consult Griffin [2] for the 
Semuvasion of this result. 

There was no data available to determine the precise 
manner in which OPARS jobs enter the input queue and so a 
Poisson arrival rate is as plausible as any other. However 
the service distribution (the interdeparture times from the 
input queue) can be measured by observing a real time display 
Of HAL's input queue and timing departures. The resulting 
data (Figure 8) and histogram (Figure 9) are surprisingly 
close to the exponential distribution. In fact, the Chi- 
Squared goodnessof fit test resulted in a test statistic of 
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Figure 8. Input Queue Interdeparture Time 
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megeo, Which indicates that the data is very close to being 
exponential with do =a. de per Nour. 

Mine final consideration is the priority ranking. 
As stated earlier an OPARS job is assigned a priority of 60 
while FNOC jobs can have priorities from 1 to 77. However, 
the batch jobs of lower priority can be ignored and we are 
Pepeewron only higher priority jobs to consider. Available 
data (Figure 10) showed higher priority jobs arriving at the 
rate of 3.45 per hour. The one remaining assumption is that 
these higher priority jobs also arrive according to a 
Bewsson Process. 


The expression for OPARS expected input queue is 


then, 
AePUT QUEUE) ] = —— —__+—__ (2) 
el Ao + AB ( 1 = Ag) 


mi 


where wu is the service rate, u is the OPARS arrival rate and 
ha memome Nhigher priority job arrival rate. This expression 
is valid only when: 


See => sd. 


u 
This input queue waiting time was computed for OPARS arrival 


rates of 2, 4, 6,...,20 per hour with results listed in 


column three of Figure 16. 


3. The Execution Sub-sSystems 


In the third and final section, the execution 
sub-system, we will characterize the OPARS program run time as 


well as delays due to other factors. 
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Because of the central memory requirements of the OPARS 
program (150K) only two OPARS jobs are allowed in the computer 
at one time. Any additional requests arriving during such 
a state would be required to wait in the input queue and would 
be passed over by the scheduling routine as it selected jobs 
for processing. Because the two OPARS jobs can be processed 
simultaneously, a queueing system of tne form GI/G/2 is 
suggested. This type of system is uncomfortable to work 
With, so simplifying assumptions will be made. 

Morse, because Of a lack of data to the contrary, 
Weewitt again ssume Poisson arrivals of OPARS jobs. It 
now remains to describe the service distribution. Unfortunately, 
even with Poisson arrivals, multi-server systems with any but 
expmential service times, are difficult to analyze math- 
ematically. 

Figures 11 and 12 contain actual OPARS run-time data 
for Low Demand and High Demand periods respectively. These 
data reflect the time an OPARS job spends in the system 
Peometo’s transfer to the execute queue until it's completion. 
Histograms of these data are in Figures 13 and 14. and are 
clearly not exponential . However, if the data can be 
meg@enaered to be from an Erlang distribution, then the nomo- 
graph in Figure 15 (see Hillier and Leiberman [4]) can be 
used to obtain an approximate value for N, the mean number 
Of OPARS jobs in the system. From this the mean service time 
can be computed via Little's result: 


mieCron |) =~ (3) 


0 
Pieeeex, +s, as before, the arrival rate for OPARS jobs. 
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Heeure i>. Nomegraph of the Long Run Average 


Number of Customers, N, in an M/E_,/2 Queue 
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A,  &=[Tsetup)] El[TCinput) J p/N/E[T(run)] —-E [T(services)] 
2/HR Souk 1.33 ny 227 500 12.43 
4/HR ee 1.39 eV 6-0 12.49 
6/HR 5.1 1.45 .3/.64/6.4 12.95 
8/HR 5.1 nee oye 13.72 
10/HR Bee 1.60 Boy ales fic 14.5 
12/HR 5.1 1.68 nent /845 15.28 
14/HR ee ae Ke “UP Vales gee 
16/HR Bal 1.88 ey ey/iake 21.28 
18/HR De. 1.99 PO /s0/ 2523 Be.359 
20/HR ee. eG 1.0/o / @ eo 


Figure 16. Busy System Summary (Time in Minutes) 


r E[T(setup) ] E[T(input queue)] p /N/ E[T(mm)] EfT(Services)] 


7 


2/HR Be 0 .O74/.15/4.5 9.6 
4/HR Biel 0 ~148/.30/4.5 9.6 
6/HR Deal 0 ~222/.46/4.6 oi 
8/HR ee 0 .296/.63/4.7 9.8 
10/HR a 0 . 369/.83/5.0 pea 
27H De 0 ~443/1.1/5.5 10/6 
14/HR Ded 0 pail el 0.0 ide 
MoyvHR 85.1 0 -591/1.7/6.4 le) 
18/HR it 0 Boe 2.0 V2 
20/HR ek 0 .738/2.8/8.4 15 


Figure 17. Idle System Summary (Time in Minutes) 





The run data was parameterized as an Erlang 


distribution with the following results: 


BUSY IDLE | 
r 4 5 
r . 6645 Rowe 
i Aenea 4.43 min 
Be CO 


Figure 18. Erlang Distribution Parameters 
Where r fanoie mean rum time. 

Mone livessthe expected execution service time, 
E{T(runs)], was computed for a series of OPARS arrivals, 
as before, with the results in column 4 of Figures 16 and 17. 

4. Summary 

In summary we have now succeeded in representing 
the total OPARS service time as the sum of three, readily 
computable sub-system times (Figure 19): 
mi@eservyice)] = E[T(setup)] + ELT(input queue)] +tE[T(run) ]. (4) 
Erlang's formula can now be employed to compute expected 


state probabilities for lines busy. 
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eee lowe SOTMULATION 


As an alternative solution technique, two computer 
programs were written to simulate the OPARS system, one each 
for High and Low Demand computer states. Output of the 
simulations were, as in the analytic model, state probabilities 
for lines busy under various demand states. 

There are two key differences in the logic flow for the 
BWomsamulations. First OPARS requests arriving during the 
busy state are sent to the input queue where they must wait 
momeemocrved. Under idle conditions jobs go directly to 
M@emeemouLer. oecondly, in the busy state, higher priority 
NObS are being created which interfere with the processing 
emer Joos. The following is a list of the various 
distributions and their parametric values used in the 


Siielay ion. 


OPARS Interarrival time: Gwe eto 4. 20/HR) 
High Priority Jobs Interarrival time: . EAP(),= 3.45/HR) 
IRG service times: NORMAL (v= 5.1 Min., o= 1) 

Input queue Interdeparture: times: =XP(\= 53.5/HR) 

OPARS run times; Idle: Govier y= (eee 34) 

OPARS run times., Busy: GAMMA (A= .701, r = 3.1) 


Beemre 20: Probability Distribution Used in Simulation 
Program 


The simulations were run for 4& hours for each OPARS 


demand rate and under each computer state. 
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meemeonly Significant factor included in the simulations 
mueriov in the analytic model had to do with a retrial 
population. In the event of all lines busy, the analytic 
model assumes arriving requests disappear, but Tieiea Lacy. 
people, when faced with a busy signal, hang up and try 
again. The simulation retries all balked attempts after 


ten minutes. 





V. RESULTS AND ANALYSIS 


Erlang formula (1) state probabilities, pi, were 
computed for OPARS request rates of 2,4,6,...,20 per hour 
under High and Low Demand conditions. Appendix A contains 
the FORTRAN program used to calculate these state probabilities 
given demand rate and mean service time, E[T(service) ]. 
Similarly Appendices D and C contain the SIMSCRIPT programs 
which simulate the OPARS system. 

ihe results of both methods are in Figures 2. and 22. 
The state probabilities are suprisingly close (differing 
by less than 2% in most cases) and clearly show that with 
eleven operating telephone lines, demand rates must be far 
in excess of projected requirements (Section III.C), before 
Deemer obabpility of all lines being busy is of concern. 

One should be aware, however, that in spite of the 
closeness of the results these two solutions merely serve 
to verify each other. Because the OPARS system is not fully 
mimeotoral at the time of this report, there is no way 
to validate either method against the actual system. 

iieaddition, if any significant changes are made to 
either the OPARS system or to the FNOC computer system, 


a new evaluation would be required. 
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APPENDIX A 


BereelTER PROGRAM TO CALCULATE ERLANG FORMULA STATE PROBABILITIES 


20 


65 
70 


90 


DIMENSION PROB (12,10), Serve (10) 
ns O001 

NN = 1 

mee, 10) (SERVE(I),I = 1,10) 
momar (rF10.5) 


mie = 12 
fem J = 1,10 

RHO = 2.0 * J/SERVE(J) 
DENOM = 0.0 

pemso K = 1,LINE 

oom = 1.0 


Memeo Lb = 1,K 
PACT = FACT * (L-1) 

Maen EQ.1) FACT = 1.0 

CONTINUE 

Dem>0 [ = 1,LINE 

FACTNU = 1.0 

poms L = 1.1 

maermyu = FACTNU * (L-1) 

miemeenG. %) FACTNU = 1.0 

CONTINUE 

PROB (I,J) = (RHO**(I-1)/PACTNU) /DENOM 
mommernOs(l,J) .LT. ZL) PROB(I,J) = 0.0 
CONTINUE 

CONTINUE 

WRITE (6,90) 

pem7o 1 = 1, LINE 

M = I-l 

WRITE(6,65) M, (PROB(I,J), J = 1,10) 
FORMAT (2X,12,10F11.6) 

CONTINUE 

IF (NN.EQ.2) STOP 

NN = 2 

GO TO 5 

FORMAT ('1) 

END 
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APPENDIX B 


SIMULATION PROGRAM: HIGH DEMAND STATE 


PREAMBLE 


MAIN 


END 


EVENT 


END 


EVENT 


END 


NORMALLY MODE IS INTEGER 

EVENT NOTICES INCLUDE HIGH.PRI.JOB, OPARS.JOB, RETRY, 
INPUT. QUEUE, 

TEMPORARY ENTITIES 

EVERY JOB HAS A TIME.OF.ENTRY AND A PRIORITY AND MAY 
BELONG TO THE QUEUE 

DEFINE QUEUE AS A FIFO SET RANKED BY PRIORITY 

THE SYSTEM OWNS THE QUEUE 

DEFINE IRG, THRUPUT, QUEUE.TIME AND TIME.OF ENTRY AS 
REAL VARIABLES 

DEFINE NO.REQUESTS, NO.ATTEMPTS, PRIORITY, LINES. 

BUSY, RUN, NO.BATCH, SYSTEM, BUSY AND EMPTY AS VARIABLES 
DEFINE SUMMARY AS A 2-DIMENSIONAL, REAL ARRAY 
ACCUMULATE STATE.PROB(O to 11 BY 1) AS THE HISTOGRAM OF 
EINES ,BUSY ~ 

ACCUMULATE QUEUE.STATS(0O TO 11 BY 1) AS THE HISTOGRAM 
OF N.QUEUE 

TALLY AVE.QUEUE.TIME AS THE AVERAGE OF QUEUE. TIME 


RESERVE SUMMARY(*,*) AS 12 BY 10 

LET RUN = 1 

LET BUSY = 1 

SCHEDULE AN OPARS.JOB NOW 

BeeEDULE A STATISTICS IN 2865 MINUTES 
START SIMULATION 


OPARS.JOB 
SCHEDULE AN OPARS.JOB IN EXPONENTIAL.F(30./RUN,4) MINUTES 
ADD 1 TO NO.ATTEMPTS 

IF LINES.BUSY EQUALS 11 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

hopper TO LINES. BUSY 

ADD 1 TO NO.REQUESTS 

LET IRG = NORMAL.F(5.1,1.0,5) 

SCHEDULE AN INPUT.QUEUE IN IRG MINUTES 

RETURN 


HEGHiotmd. JCB 

SCHEDULE A HIGH.PRI.JOB IN EXPONENTIAL.F(17.4,3) MINUTES 
CREATE A JOB 

EPP RuORITY (JOB) = 2 

FILE JOB IN QUEUE 

RETURN 
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EVENT 


END 


EVENT: 


END 


EVENT 


END 


RETRY 
IF LINES. BUSY EQUALS 11 

ADD 1 to NO.ATTEMPTS 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO NO.REQUEST 

Aep) i TO LINES. BUSY 

LET IRG = NORMA.F(5.1,1.0,5) 

SCHEDULE AN INPUT.QUEUE INIRG MINUTES 
RETURN 


INPUT. QUEUE 
IF SYSTEM EQUALS EMPTY AND NO.BATCH IS LESS THAN 2 
SCHEDULE A JOB.COMPLETION IN GAMMA.F(4.43,3.1,7)MINUTES 
ADD 1 TO NO.BATCH 

OTHERWISE 

CREATE A JOB 

LET TIME.OF.ENTRY(JOB) = TIME.V 

fee PRLORITY(JOB) = 1 

FILE JOB IN QUEUE 

REGARDLESS 

RETURN 


EXECUTION 

IF SYSTEM EQUALS BUSY 

SCHEDULE AN EXECUTUION IN EXPONENTIAL.F(1.12,2) MINUTES 
REGARDLESS 

IF N.QUEUE EQUALS O 

RETURN 

OTHERWISE 

i PREORITY(F.QUEUE) = 2 

REMOVE FIRST JOB FROM QUEUE 

DESTROY THE JOB 

RETURN 

OTHERWISE 

IF NO.BATCH . EQUALS 2 

RETURN 

OTHERWISE 

REMOVE FIRST JOB FROM QUEUE 

LET QUEUE.TIME= (TIME.V -TIME.OF.ENTRY(JOB) * 1440. 
DESTROY THE JOB 

fier oraM EQUALS BUSY 

SCHEDULE A JOB.COMPLETION IN GAMMA.F(6.06,4.34,6)MINUTES 
OTHERWISE 

SCHEDULE A JOB.COMPLETION IN GAMMA.F(4.43,3.1,7)MINUTES 
REGARDLESS 

ADD 1 TO NO.BATCH 

RETURN 
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EVENT JOB.COMPLETION 
DUBTRACT 1 FROM LINES.BUSY 
SUBTRACT 1 FROM NO.BATCH 
IF N.QUEUE EQUALS 0 
RETURN 
OTHERWISE 
IF SYSTEM EQUALS BUSY 
CANCEL THE EXECUTION 
REGARDLESS 
SCHEDULE AN EXECUTION NOW 
RETURN 

END 


EVENT STATISTICS 
START NEW PAGE 
PRINT 1 LINE WITH RUN*2 THUS 
ARRIVAL RATE = **HOUR 
SKIP 3 OUTPUT LINES 
PRINT 1 LINE THUS 
PR(BUSY LINES) PR(INPUT QUEUE) 
SKIP 2 OUTPUT LINES 
Boner = 1 TO 12, DO 
PRINT 1 LINE WITH I-1, STATE.PROB(I)/2 AND QUEUE. 
pooms(l)/2 THUS 
¥¥ KHEXHE SEXES 
LOOP 
SKIP 5 OUTPUT LINES 
PRINT 3 LINES WITH AVE.QUEUE.TIME,NO.ATTEMPTS AND NO. 
REQUESTS THUS 
AVERAGE OPARS QUEUE TIME = #%#, ###* 
Veom@osrs . ATTEMPTED RRRE 
NO.OPARS COMPLETE *#%#% 
Omer = 1 to 12, DO 
LET SUMMARY(I,RUN) = STATE.PROB(I)/2 
LOOP 
SCHEDULE A STATISTICS IN 2880 MINUTES 
RESET TOTALS OF LINES.BUSY, N.QUEUE AND QUEUE.TIME 
IF RUN EQUALS 10 
GO TO FINAL.STATS 
OTHERWISE 
ADD 1 TO RUN 
LET NO.ATTEMPTS ° 
PeleNoO. REQUESTS 
RETURN 
'PINAL.STATS' 
START NEW PAGE 
PRINT 5 LINES THUS 
EMPTY SYSTEM SUMMARY 
STATE 2/HR 4/HR 6/HR 8/HR 10/HR 12/HR 14/HR 16/HR 18/HR 20/HR 
Henwr = 1 to 12, DO 
PRINT 1 LINE WITH I-1, SUMMARY (I,1), SUMMARY (1,2), 
SUMMARY(I,3), SUMMARY (1,4),SUMMARY(1I,5),SUMMARY(I,6), 
SUMMARY(1I,7), SUMMARY (1,8) ,SUMMARY(I,9),SUMMARY(I,10) ,THUS 
¥*¥ , RREE , ERRE , BREE , KEKE , RRER we RNE , RRRE , EHRE _RRKRE” HERE 
LOOP 
STOP 
END 46 
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APPENDIX C 


SIMULATION PROGRAM: LOW DEMAND 


PREAMBLE 


END 


MAIN 


END 


EVENT 


END 
EVENT 


NORMALLY MODE IS INTEGER 

foe NOTICES PNCLUDE HIGH.PRT.JOB,OPARS.JOB,RETRY,INPUT. 
QUEUE, 

ec ORARY ENTITIES 

BvenhyY JOB HAS A TIME.OF ENTRY AND A PRIORITY AND MAY 
BELONG TO THE QUEUE 

BEPINE QUEUE AS A FIFO SET RANKED BY PRIORITY 

THE SYSTEM OWNS THE QUEUE 

Peek LRG, THRUPUT,QUEUE.TIMs AND TIME.OF.ENTRY As 

REAL VARIABLES 

Pee NOoREQUBSTS; NO. ATTEMPTS ,PRIORITY, LINES.BUSY, 
RUN, NO.BATCH, SYSTEM, BUSY AND EMPTY AS VARIABLES 
DEFINE SUMMARY AS A 2-DIMENSIONAL, REAL ARRAY 
PeeOMULATE STATE.PROB(O TO 11 BY 1) AS THE HISTOGRAM 

Si otNES.BUSY 

PeouMUBATE QUEUE.STATS.(O TO 11 BY 1) AS THE HISTOGRAM 
OF N.QUEUE TALLY AVE.QUEUE.TIME AS THE AVERAGE OF QUEUE. 
PME 


RESERVE SUMMARY(*,*) AS 12 BY 10 

eT RUN = 1 

LET SYSTEM = 1 

mee BUSY = 1 

SCHEDULE AN OPARS.JOB NOW 

SCHEDULE AN EXECUTION IN EXPONENTIAL.F(1,12,2)MINUTES 
SCHEDULE A HIGH.PRI.JOB IN EXPONENTIAL.F(17.4,3)MINUTES 
CREATE A JOB 

item PREORITY (JOB) = 2 

FILE JOB IN QUEUE 

SCHEDULE A STATISTICS IN 2880 MINUTES 

START SIMULATION 


OPARS.JOB 

SCHEDULE AND OPARS.JOB IN EXPONENTIAL.F(30./RUN,4) MINUTES 
meee TO NO. ATTEMPTS 

IF LINES.BUSY EQUALS 11 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

Peper tO LINES.BUSY 

PDO NO. REQUESTS 

metewRe = NORMA.F(5.1,1.0,5) 

SCHEDULE AN INPUT.QUEUE IN IRG MINUTES 
RETURN 


HIGH.PRI.JOB 
SCHEDULE A HIG.PRI.JOB IN EXPONENTIAL.F(17.4,3)MINUTES 


CREATE A JOB 
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APPENDIX C 


SIMULATION PROGRAM: LOW DEMAND STATE CONT'D 


END 


LET PRIORITY(JOB) = 2 
FILE JOB IN QUEUE 
RETURN 


EVENT RETRY 


END 


IF LINES.BUSY EQUALS 11 

ADD 1 TO NO.ATTEMPTS 

SCHEDULE A RETRY IN 10 MINUTES 

RETURN 

OTHERWISE 

ADD 1 TO NO.REQUESTS 

hep 1 TO LINES.BUSY 

mere tRG = NORMAL.F(5.1,1.0,5) 

SCHEDULE AN INPUT.QUEUE IN IRG MINUTES 
RETURN 


EVENT INPUT. QUEUE 


END 


IF SYSTEM EQUALS EMPTY AND NC.BATCH IS LESS THAN 2 
SCHEDULE A JOB.COMPLETION IN GAMMA.F(4.43,3.1,7) MINUTES 
feo 2 TO NO.BATCH 

OTHERWISE 

CREATE A JOB 

LET TIME.OF.ENTRY(JOB) = TIME.V 

for PRIORITY(JOB) = 1 

FILE JOB IN QUEUE 

REGARDLESS 

RETURN 


EVENT EXECUTION 


mor oTEM EQUALS BUSY 

SCHEDULE AN EXECUTION IN EXPONENTIAL.F(1.12,2) MINUTES 
REGARDLESS 

iN. QUEUB EQUALS C 

RETURN 

OTHERWISE 

Eee REORITY (F.QUEUE) = 2 

REMOVE FIRST JOB FROM QUEUE 

DESTORY THE JOB 

RETURN 

OTHERWISE 

IF NO.BATCH EQUALS 2 

RETURN 

OTHERWISE 

REMOVE FIRST JOB FROM QUEUE 

LET QUEUE.TIME = (TIME.V - TIME.OF.ENTRY(JOB)) * 1440. 
DESTORY THE JOB 

Peo c STEM EQUALS BUSY 

SCHEDULE A JOB.COMPLETION IN GAMMA.F(6.06,4.34,6)MINUTES 
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SIMULATION PROGRAM: LOW DEMAND STATE CONT'D 


OTHERWISE 

SCHEDULE A JOB.COMPLETION IN GAMMA.F(4.43,3.1,7) MINUTES 
REGARDLESS 

ADD 1 TO NO.BATCH 

RETURN 


END 


EVENT JOB.COMPLETION 
SUBTRACT 1 FROM LINES .BUSY 
SUBTRACT 1 FROM NO.BATCH 
If N.QUEVE EQUALS 0 
RETURN 
OTHERWISE 
ZieotolEM EQUALS BUSY 
CANCEL THE EXECUTION 
REGARDLESS 
SCHEDULE AND EXECUTION NOW 
RETURN 

END 


Se rewoLATISTICS 
START NEW PAGE 
Boe | LENE WITH RUN*e THUS 
ARRIVAL RATE = **/HOUR 
mere 3 OUTPUT LINES 
eaeNe | LINE THUS 
PR(BUSY LINES) PR(INPUT LINES) 
mere Cc OUTPUT LINES 
monet = 1 TO l2, DO 
Reese EINE WITH IT-1, STATE.PROBLE(I).2 AND QUEUE.STATS 
@i)72 THUS 
x * REESE RRR 
[EGOP 
eee > 6 OUTPUT LINES 
Beene 35 LINES WITH AVE.QUEUE.TIME.NO.ATTEMPTS AND NO. 
BeQunots THUS 
AVERAGEOPARS QUEUE TIME = #%, ###% 
NO.OPARS ATTEMPTED *#%#%# 
NO OPARS COMPLETED ¥**%*%# 
mon ©: = L TO le, DO 
fet SUMMARY(IL,RUN) = STATE.PROB(I)/2 
LOOP 
SCHEDULE A STATISTICS IN 2880 MINUTES 
fee hOrnio OF LINES.BUSY,.N.QUEUE AND QUEUE.TIME 
ma RUN EQUALS 10 
come PENAL.STATS 
Sine RWISE 
mom 1 TO RUN 
Per NO.ATTEMPTS = 0 
fer NO. REQUESTS = Q 
mi URN 
SeENAL.oLATS ' 
START NEW PAGE 
Pees LENES THUS 
Sus. of STEM SUMMARY 
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Peete /7HR 4/HR 6/HR 8/HR 10/HR 12/HR 14/HR 16/HR 18/HR 20/HR 


HOR It = 1 TO 12, DC 

PRINT 1 LINE WITH I-1, SUMMARY (I,1), SUMMARY (I,2), 
SUMMARY (1,3) ,SUMMARY(1I,4),SUMMARY(I,5) ,SUMMARY(I,6), 
SUMMARY (1,7), SUMMARY(I,8), SUMMARY (1,9), SUMMARY 


10) THUS 
RH  KRKKK  KHKK  KHKK  RKKK KKK KKK KKK  KKKK  KHKK HHH 


LOOP 
STOP 
END 
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